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ABSTRACT 



An experimental investigation was conducted to examine 
convection heat transfer in a duct of rectangular cross-sec- 
tion, having both a straight and curved flow passage. One 
wall was kept at a constant temperature, with the opposite 
wall being adiabatic. Heating of the air flowing in the 
channel was accomplished through Joulean heating of a wall 
composed of Temsheet. 

The experiments were conducted for steady state in 
Laminar and Transition flows. Development and progagation of 
the Taylor-Gortler vortices was shown to enhance the heat 
transfer rate in the curved section as compared to the 
straight section. 
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I. INTRODUCTION 



A. TAYLOR-GORTLER VORTICES 

Since the early part of this century, considerable data 
and numerous studies have demonstrated that fully developed 
laminar flow along a concave wall does not remain two-dimen- 
sional [Refs. 1, 2, 3J. Instead, the flow forms a system of 
spiral vortices, of counter rotating pairs, whose axes are 
aligned in the direction of the principle fluid flow. This 
phenomenon, known as Taylor-Gortler vortices, is the result 
of the variations in the centifugal forces acting on the 
fluid particles. Schematics of this type of motion are shown 
in Figures 1 and 2. 

In a channel that is curved in the streamwise direction, 
the fluid particles near the center of the flow cross-section 
are subjected to higher centifugal forces than those slower 
moving particles nearer to the boundary wall. Consequently, 
the tendency is for the fluid in or near the center of the 
channel to move outwardly towards the concave wall. This 
action then causes the fluid particles near the boundary wall 
to move in the spanwise direction, and finally move radially 
inward, replacing the outwardly moving particles. Once in the 
center of the channel these fluid particles come under those 
same stronger forces of the free stream velocity and continue 
the rotational effect. This cyclic motion' causes the 
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Figure 1. Schematic of Taylor vortices between cylinders 
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Primary Flow Direction 




Figure 2. Schematic of Taylor-Gortler vortices between 
parallel plates 
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formation, and propagation, of th.e counter-rotating Taylor- 
Gortler vortices. 

Although the Taylor-Gortler vortices are primarily a 
laminar flow phenomenon they may affect the transition from 
laminar to turbulent flows, and are therefore important to 
understand {Ref. 2]. The additional mixing provided by the 
secondary motion has been thought to account for the increase 
in the rate of heat transfer from concave curved walls versus 
straight walls [Ref. 4 ] . 

The Taylor-Gortler vortices appear to have similarities 
to other vortex flow patterns such as the longitudinal vor- 
tex rolls developed in the laminar forced convection heating 
of fluid layers between parallel plates [Ref. 5], It has 
also been thought that the cross-hatching, observed in 
reentry vehicles, can be explained at least in part, by the 
presence of streamwise vortices [Ref. 6] . 

There are many possible applications that could result 
from a more thorough understanding of the Taylor-Gortler 
vortices and their effect on heat transfer and fluid flow 
characteristics. Applications, such as improved turbine 
blade cooling [Refs. 7, 8, 9], heat exchanger designs, as well 
as other such engineering designs and capabilities could take 
advantage of these vortices and improved heat transfer 
characteristics of the flow in the presence of these vortices. 
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B. HISTORY 



The instability of an inviscid fluid rotating symmetri- 
cally about an axis, as in curved flow, was first considered 
by Lord Rayleigh in 1916 [Ref. 10]. In assuming that the 
fluid was non-viscous, he determined that a necessary 
condition, for the stability of an inviscid fluid motion, was 
that the circulation increased with increasing radius. G. I. 
Taylor, in 1923, [Refs. 1, 11] expanded the earlier works of 
Rayleigh by an extensive analytical and experimental study of 
viscous fluids. His investigations focused on the flow 
between two cylinders, in which the inner cylinder rotated 
while the outer cylinder remained stationary. Taylor 
ascertained that such couette flows became unstable when the 
value of the dimensionless Taylor number, defined as: 




exceeded a critical value of 41.3. In the above equation, 
*d' is the width of the gap, assumed small when compared to 
' R^ ' , the radius of the inner cylinder, 'U^' is the periphal 
velocity of the inner cylinder, and 'Re' is the Reynolds 
number. Taylor determined that for those cases in which the 
value of the Taylor number became greater than the critical 
value, a secondary motion developes and the Taylor vortices 
form. 
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Instability of a similiar nature is also observed when a 
viscous fluid flows in a curved channel due to a pressure 
gradient acting along the channel wall. This problem was 
considered analytically first by W. R. Dean [Ref. 12J in 
1928, for a channel formed by two concentric cylinders, 
where the radius of the inner cylinder was large in comparison 
to the small spacing between the inner and outer cylinder 
walls. Dean, in his studies, concluded that there would be 
an initiation of the flow instability, and the formation of 
vortices, (similiar to the Taylor vortices), when the Dean 
number, defined as: 

De = Re 

achieved a value greater than 36. In the above equation, the 
Reynolds number is based on the mean velocity of the unper- 
turbed flow, 'd* represents the channel half-width, and ' R^ ' 
is again the inner cylinder radius. The analytical work of 
Dean, was later verified by W. H. Reid [Ref. 13] , using an 
approximate numerical solution. 

H. Gortler [Ref. 2], in 1940, studied the stability of 
laminar boundary layer profiles, on curved walls under the 
influence of small disturbances. In his study, Gortler found 
that these disturbances were similiar in nature to those 
studied by G. I. Taylor, which had led to the fluid flow 
instability in the form of vortices. By approximate numerical 
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calculations , Gortler concluded that only on the concave 
boundary walls were th.e amplified disturbances produced, and 
that the overall flow profile appeared to remain laminar in 
nature. In 1955, Gortler's approximate results were verified, 
with an exact solution, by G. Hammerlin, as reported by H. 
Schlichting [Ref. 14] and then further substantiated in a 
more extensive numerical analysis by A. M. 0. Smith [Ref. 3]. 
More recently, verification of these numerical solutions have 
been demonstrated by use of hot wire anemometry, laser doppler 
systems, and flow visualization techniques, [Refs. 15, 16], 
and in 1976, Y. Aihara [Ref. 17] conducted a non-linear 
analysis of the Gortler vortices. 

With the growing interest in the effects of the secondary 
flows associated with the Tay lor- Gor tl er vortices, numerous 
specific studies have been published concerning the influences 
of these vortices on the transfer of heat in both the laminar 
and turbulent flow regimes. F. Kreith [Ref. 4], in 1955, 
studied the influence of heat transfer with respect to the 
curvature of the boundary wall for fully turbulent flows, 
and concluded that the heat transfer from the heated concave 
boundary wall was considerably higher than that transferred 
from the convex boundary wall of the same curvature and under 
similiar turbulent flow conditions. 

In 1965, L. Persen [Ref. 18], considering the special 
cases of very high and very low Prandtl number fluids , 
related the increase in heat transfer rate from a curved wall 
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to the presence of the Taylor- Gortler yortices. His con- 
clusions were that the effect, of the development of the 
vortices, was to increase the rate of heat transfer through 
the boundary layer. 

One of the first experimental works on the effects of the 
Taylor-Gortler vortices with respect to enhanced heat transfer 
was by P. McCormack, et al., [Ref. 19], published in 1970. 

With the data, they concluded that it was necessary to retain 
the non-linear terms of the flow equation in order that the 
complete theoretical solution of the effects of the Taylor- 
Gortler vortices on heat transfer be explained. In 1977, 

R. Kahawita and R. Meroney [Ref. 20], furthered McCormack’s 
work and concluded that the higher order terms, as well as 
the normal velocity components of the main flow, became 
increasingly more significant to the calculations at smaller 
wave numbers. 

In considering the forced convective heat transfer for 
fully turbulent flows, between horizonal parallel flat plates, 
Y. Mori and Y. Uchida [Ref. 5], in 1966, reported the 
formation of longitudinal vortex rolls with an axis parallel 
to the flow direction. These streamwise-axis directional 
rolls, similiar to the Taylor-Gortler vortices, formed when 
a critical value of temperature difference between the 
plates was exceeded. In 1971, M. Akiyama, et al., [Ref. 21], 
confirmed the formation of these vortices at a critical 
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Rayleigh number of 17Q8, for the case of heating from the 
bottom plate surface. 

There has been a rather limited amount of published 
literature regarding the flow and heat transfer in curved 
channnels of rectangular cross-section. Much of the pub- 
lished literature involves the development of numerical 
approximations and solutions for heat and mass transfer in 
curved ducts of various geometries (most notably, that of 
circular and triangular ducts) . K. Cheng and M. Akiyama 
[Ref. 22] developed a numerical solution for forced 
convection heat transfer with laminar flows in curved channels 
of rectangular cross-section, but only for small aspect 
ratios, (B/D; where, 'B* is the channel width and ’ D * is the 
depth of fluid flow) . Using the fully elliptical forms of 
the transport equations, G. Yee and J. A. C. Humphrey 
[Ref. 23] , developed a numerical solution for flows and 
boundary conditions similiar to Cheng and Akiyama. K. Cheng, 
et al., [Ref. 24], in 1974, using numerical methods solved the 
Graetz problem for a curved channel of square cross-section 
CB/D=1) . 

In 1963, B. S. Petukhov and V. N. Popov [Ref. 25], obtained 
an analytical expression for the Nusselt number for the flow 
of an incompressible fluid with variable properties in a 
circular tube. In 1976, A. A. Shibani and M. N. Ozisik 
[Ref. 26] using matched asymptotic expansion techniques for 
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the case of uniform wall temperature, solyed the heat 
transfer problem between parallel plates with turbulent flow, 
for a wide range of Prandtl numbers. 

As to analytical and experimental studies of the parallel 
plate, channel flow heat transfer problem, Y. Mori, et al., 
[Ref. 27] , obtained results for hydrodynamically fully 
developed flows with constant wall heat flux, in curved 
channels of square cross-section. W. M. Kays and E. Y. 

Leung, [Ref. 28], experimentally obtained solutions for 
turbulent flow heat transfer in a concentric circular tube 
annulus with fully developed velocity profile and constant 
heat rate per unit length, for a fluid of Prandtl number 0.7. 
M. Durao [Ref. 29] and J. Ballard [Ref. 30], reported results 
for large aspect ratio (B/D = 40) rectangular cross-sectional 
channels with laminar flows. P. F. Brinich and R. W. Graham 
[Ref. 31] , reported results for turbulent flows in a rectang- 
ular curved channel with an aspect ratio of 6, for the inner 
wall heated, the outer wall heated, and both walls heated. 

In 1979, F. D. Haynes and G. D. Ashton [Ref. 32], in studying 
heat transfer from a river to its ice cover, reported results 
for rectangular cross-section channels of aspect ratio 10, 
and compared them to the data obtained by Ashton [Ref. 33] , 
and Hsu [Ref. 34] . 
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IT. INTENT OF THE STUDY 



The primary purpose of this study was to investigate the 
effect of the Taylor-Gortler vortices, on the enhancement of 
tfie heat transfer rate in a curved channel of rectangular 
cross-section and large aspect ratio. This investigation 
began with flow velocities that were in the laminar flow 
regime and continued into the transition region, between 
laminar and turbulent flow. Further, this study compared the 
heat ransfer rate obtained in a curved section with the rate 
obtained in a straight section of identical aspect ratio. 

It was eocpected that the presence of the Taylor-Gortler 
vortices would enhance the heat transfer process at all 
velocities in the flow regimes investigated. The cause of 
this would be the secondary flow velocity components of the 
vortices, transporting heated fluid from the concave boundary 
wall inward, while simultaneously displacing the cooler fluid 
particles in the center, causing them to move outwardly 
towards the heated concave wall. 

This study was conducted using a single rectangular cross- 
sectional channel, that incorporated both a straight and a 
curved test section. The results obtained in the straight 
section were compared to those results obtained in the curved 
section, at the same flow velocities, in an effort to 
determine the effects of the Taylor-Gortler vortices on the 
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transfer of heat. Additionally, the results of both the 
straight and curved section tests were compared to the 
experimental results of Ballard jRef. 30] and Durao [Ref. 29], 
Those results of the straight section that fell into the 
laminar flow region, were also compared to the analytical 
solution as given by P. McCuen, et al . , [Ref. 35]; while 
those in the transition region were compared to the 
experimental results of Kays and Lueng [Ref. 28] , and the 
analytical results of the Dittus -Boelter equation [Ref. 38] 
and Eckert's equation [Ref. 39]. While these equations are 
for fully turbulent flows in circular tubes they have been 
used to compare turbulent flows in other geometry channels 
when the diameter of the tube is replaced by the hydraulic 
diameter (wetted permimeter divided by the cross-sectional 
area) . 

The straight section results, for transition flow, were 
also compared to the analytical solutions of Petukhov and 
Popov [Ref. 25] and Shibani and Ozisik [Ref. 26], the former 
for circular tubes while the latter is for parallel plates, 
both with fully turbulent flows. 

The curved section results for flows in the transition 
region were compared to the experimental results of Brinich 
and Graham [Ref. 31], and also to the solutions of Petukhov 
and Popov, and of Shibani and Ozisik, mentioned above. 
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III. EXPERIMENTAL WORK 



A. DESCRIPTION OF THE APPARATUS 

A channel of rectangular cross-section built as described 
in references 29 and 30, was used to achieve the objectives 
of this study. The channel was manufactured from two 0.635 
centimeter thick sheets of plexiglas, separated by 0.635 
centimeter spacers whose inside lengths were used as the 
sides of the channel. 

The channel, as shown in Figure 3, was composed of a 
straight section, 122.0 centimeters long, followed by a 
curved section of 180 degree arc. The interior concave wall 
of the curved section has a radius of curvature of 30.5 
centimeters. The rectangular channel, with an aspect ratio 
of 40, was 0.635 centimeters high by 25.4 centimeters wide, 
with a cross-sectional area of 16.15 square centimeters and 
a wetted perimeter of 57.07 centimeters. A cross-sectional 
view of the channel, for both straight and curved sections, 
is shown in Figure 4. 

The working fluid used in this investigation was air, at 
room temperature, which entered the channel via an entrance 
bell constructed of Plexiglas, which was connected to the 
straight section. The entrance bell was designed and 
manufactured in accordance with ASME nozzle standards, with 
an elliptical curved base on a major axis equal to ten inches 
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Inlet Nozzle Straight Test Section 
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Plexiglas 

Aluminum 




Figure 4. Cross-Sectional view of Channel 
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and a minor axis of one inch. In an attempt to prevent the 
introduction of dust and other foreign matter into the flow 
channel, cheese-cloth was attached to the opening of the 
entrance bell nozzle. After passing through the channel, 
the flow was directed through an aluminum exhaust nozzle 
which was connected, by 3/4 inch PVC piping, to a Fisher and 
Porter Company variable area flow meter model 10A3565A. 

The rotameter, with a Bead Guide Tube number FP-1-35-G-10 
and a 1GSVGT float serial number 68T60, at 100% full scale 
had a flow rate of 0.793 cubic meters of air per minute 
(28.0 standard cubic feet of air per minute at one atmosphere 
and 70 0 F) . 

The flow of air was drawn through the channel, piping and 
rotameter by a single electrically driven Spencer Turbo 
Compressor, 30 Hp. at 3500 rpm, rated at 550 cubic feet per 
minute at 70°F and one atmosphere. The flow rate was con- 
trolled by the use of valving both at the compressor and at 
the exit of the rotameter. 

To obtain the experimental heat transfer data for this 
study, two test sections of the channel were constructed: 

(.1) a straight test section, 29.2 centimeters in length, 
located far enough downstream to ensure hydrodynamically 
fully developed flow, with a heated test section area of 
741.9 square centimeters. (2) a curved test section, 28.3 
centimeters in length, subtending an arc of 53.1 degrees, 
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located in the lower half of the concaye curved portion of 
the channel, with a test section area of 717.4 square 
centimeters . 



In each of the test sections, the outer wall was 
modified as shown in Figures 5 and 6. Temsheet, a carbon 
impregnated porous paper with, a property of uniform 
electrical resistivity, was glued to the interior surface of 
the outer wall of each of the two test sections. The flow of 
air was then to be heated by Joulean heating of the Temsheet. 
As the electrical resistance of Temsheet is not constant, but 
rather slowly variable with temperature, a precision resistor 
with an electrical resistance of 2.078 ohms was connected in 
series with the Temsheet to allow calculation of the 
instanteous power being supplied. 

The variables to be measured, in this investigation, were 
(.1) the inlet flow temperature for each of the test sections 
(Ti N ) , (2) the outlet flow temperature at the exit of each of 
the test sections (Tq^) > (.3) the temperature of the heated 
boundary wall of each of the test sections (T ) , (4J the 
temperature of the unheated boundary wall of each of the test 
sections (.T^) , (.5) the flow temperature at the exit of the 
channel , (6) the temperature between the two insula- 

tion layers for each of the two test sections (TjxjgJ > (7) 
the voltage across the precision resistor (Vp p ) , (8) the 
heater voltage , and (9) the flow rate through the 

rotameter (Q) . 
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Figure 5. 




Straight Test Section, detailed schematic 
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FLOW 




Insulation 
Temsheet 



Plexiglas 



28.3 cm 



Thermocouples 



Figure 6. Curved Test Section, detailed schematic 
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The temperature measurements required ahoye, as well as 
all other temperature measurements in this investigation, 
were made with Copper-Constantan thermocouples. All the 
thermocouples had been previously calibrated using a 
ROSEMOUNT Communtating Bridge model 920A, and a ROSEMOUNT 
Constant Temperature Bath model 913A, as outlined by Ballard 
{Ref. 30] . 

A total of sixty-four Copper-Constantan, glass insulated, 
30 gauge thermocouples were constructed and inserted in 
specific locations to facilitate the measurement of the 
required temperatures. All the thermocouples were wired to a 
Hewlett Packard Data Acquisition System, model 2010c, which 
enabled an automatic printed record of the instantaneous 
reading for each thermocouple location. All the thermocou- 
ples were also connected to a Kaye Instrument Ice Point 
Reference, model K140-4, serial 1780. 

In each of the test sections, five of the thermocouples 
were inserted between the first and second layers of 
ARMAFLEX 22 insulation, and another five inserted between the 
second and third layers. Each of these sets of five 
thermocouples were connected in parallel to facilitate the 
recording of an average value of the temperature (TINS) 
between the respective layers of insulation. Additionally, 
sets of four thermocouples, also connected in parallel to 
indicate average temperature, were inserted into the channel 
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at the following three locations; (1) the entrance of the 
channel, (2) the exit of the straight test section, and (3) 
the exit of the curved test section. At each of these three 
locations, the thermocouples were spaced five centimeters 
apart in the spanwise direction to allow reading of the 
bulk temperature of the flowing air. 

In each of the two test sections, eight thermocouples 
were placed in direct contact with the Temsheet through 
small diameter holes drilled into the Plexiglas outer walls, 
with the thermocouple beads electrically insulated from the 
Temsheet with ENMAR Heat Resisting Glyceryl Phthalate. The 
final sixteen thermocouples were placed, eight to a test 
section in the Plexiglas unheated boundary wall, opposite 
those positioned in the Temsheet. A detailed sketch of the 
thermocouple placement is provided in Figure 7. The average 
heated outer wall temperature (T ) , was provided by those 
thermocouples in contact with the Temsheet. The average 
unheated inner wall temperature (T w ^) , was provided oy tne 
tnermocoupies embedded in the Plexiglas wall. 

Tne thermal insulation consisisted of layers of 
ARMSTRONG ARMAFLEX 22 Sheet Insulation, a flexible foamed 
plastic material, on the heated side of each test section, 
and black foam insulation one-half inches thick elsewhere. 
Three layers of the ARMAFLEX 22 insulation were used, each 
layer had a thickness of 0.636 centimeters, and covered an 
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Figure 7. Thermocouple Placement in Test Sections 
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area slightly larger than the respective test section. 

These layers of insulation were used in the compution of 
heat flux losses to the environment, sample calculations 
of these losses are given in Appendix C. The rest of the 
channel was covered with the one-half inch thick Black 
foam insulation to reduce the losses to the surroundings 
through the plexiglas walls as the heated air transited 
the channel. The insulation was held in place by the use 
of Adhesive Heat Resistant ventilation Duct Tape as seen 
in Figure 8. The entire length of the channel as well as 
all the connections between the channel and the rotameter, 
were sealed with General Electric Silicone Rubber Sealant 
Caulk to ensure that there was no leakage of air into the 
channel or rotameter downstream of the entrance bell 
nozzle . 

Two. aluminum electrodes, 0.318 centimeters thick, were 
inserted in between the Temsheet and the plexiglas, in each 
of the test sections, as shown in Figure 4. The plates, 
with the Temsheet and aluminum electrodes were assembled 
in the channel with plastic teflon screws for the purpose 
of electrical insulation. 

B. EXPERIMENTAL PROCEDURES 

Experiments were conducted with the volumeteric flow 
rate varied incrementally from 0.198 to 0.674 cubic meters 
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Figure 8. Photograph of Test Channel 
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of air per minute to correspond to particular Reynolds 
numbers. The data was recorded and calculations performed 
for each flow rate. The experimental procedures followed 
were the same for both the straight and curved test 
sections. The straight test section runs served as a 
baseline for the curved test section and allowed tor 
comparison of the results. The appropriate electrical 
power required to heat the Temsheet boundary wall to 
approximately 50°U was determined through preliminary 
experiments. This temperature was choosen to allow at 
least a twenty degree celisus temperature difference 
between the unheated and heated boundary wall average 
temperatures and to facilitate achieving steady state 
conditions within a reasonable time span. Estimates of 
time to achieve steady state conditions and decisions as 
to data collection were also based on these preliminary 
test runs. 

The determination as to what constituted a steady 
state condition required careful consideration in view 
of the multiple input variables and their effects on the 
results. The final criteria for steady state was based on 
three input variables: Cl) flow rate percentage, as 

recorded from the rotameter; (2) heated boundary wall 
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temperature, as recorded from the thermocouples in that 
area; and (3) the heater voltage, as recorded on the 
voltmeter. When these conditions varied by less that 
two percent over a ten minute period, the experiment 
was considered to have achieved steady state. Under 
these conditions the assumption of steady flow and 
constant heat flux from the heated wall were assumed 
valid . 

Originally, the time required to reach this steady 
state condition was found to be in the range of two 
hours, accordingly all runs were at two and one-half 
hours durations. Data was taken on one hour intervals 
for the first two hours and then at ten minute intervals 
for the next thirty minutes. Later it was found that the 
higher flow percentages for the curved test section, 
required about thirty minutes more time to full achieve 
steady state, so test runs with a flow percentage of 
sixty percent or greater were extended with three 
additional ten minute interval readings for a total time 
of three hours. 

The instantaneous values for all the temperatures 
were automatically printed by the digital recorder in 
millivolts. A standard thermocouple table was used in 
conjunction with previously obtained thermocouple 
calibration data to convert the millivolts to degrees 
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Celsius. The yolumetnc flow rate was read directly from 
the installed rotameter as the percent of full scale flow. 
The instantaneous power supplied was calculated using the 
relationship : 

V H V PR 

Q p = 

r pr 

Where the precision resistor voltage (Vp^) and the heater 
voltage (V^) were recorded in volts from the installed 
voltmeter in the digital Data Acquistion System, and the 
electrical resistance of the precision resistor (Rpp) had 
been previously determined to be 2„078 ohms. 

As the area of the curved test section was slightly 
less than the area of the straight test section, (0.0718 
vs. 0.0742 square meters), the heat flux from the heated 
wall, of the respective test sections was different. The 
room temperature entering the test section was checked 
with the use of a mercury thermometer located in the 
vicinity of the apparatus. 

During the preliminary test runs of the straight 
section, the thermocouple array recording the exiting 
temperature of the flow of air, was set at three 
millimeters from the heated wall. This was satisfactory 
at the lower volumetric flow rates, but at the higher 
flow rates it led to calculations in which the power 



38 



convected to the air was greater, by approximately fifteen 
to twenty percent, than the power being supplied to the 
heater. It was determined that the thermocouples were 
not indicating the correct bulk temperature of the 
exiting flow of air, but rather a higher value of temp- 
erature, this then resulted in an erroneously high value 
of convected heat transfer, and higher Nusselt numbers. 

To correct this, the system was allowed to run until the 
walls of the downstream curved section achieved steady 
state. Then, calculating the losses, by use of the 
thermocouples in between the insulation layers of the 
curved test section, and applying these losses to the 
entire curved section area, the energy loss was calculated. 
The flow exiting the channel was assumed to be thoroughly 
mixed and therefore the temperature as recorded by the 
positioned thermocouples, at the exit of the curved 
section, was used as an accurate representation of the 
average bulk exit temperature. With the final exit 
temperature and the losses now known, the temperature of 
the air as it exited the straight section was calculated. 

To correct this condition, the thermocouples were 
inserted further in to the channel, or raised closer to 
the outer wall until the millivolt reading was compatible 
with the calculated correct temperature. After this 
procedure, accomplished at a flow rate of 0.595 cubic 
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centimeters of air per minute C?5% of full scale) was 
completed, additional test runs were conducted at 0.555 
and 0.635 cubic centimeters of air per minute (701 and 
80% respectively) , without adjusting the thermocouple 
positions, to check the accuracy. In each case, the 
temperature was within 0.25°C of what the energy balance 
indicated it should be. 

It was later discovered that the curved section re- 
sults had the same problem at flow rates greater then 
0.595 cubic centimeters of air per minute. Since the 
channel construction is such that it does not allow 
further downstream data collection, and since the power 
difference was less than four percent it was decided to 
proceed with the correlation of the data and experimenta- 
tion and note the discrepancies. 
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IV. PRESENTATION OF DATA 



A. ANALYSIS 

The uniform electrical resistivity of the Temsheet, 
approximated a constant heat flux surface for the heated 
wall. With insulation on the unheated wall the heat losses 
through that wall were negligible and the wall could be 
considered adiabatic. With the large aspect ratio of 40, 
the channel was considered to be approximately parallel 
plates. The straight portion of the channel upstream of 
the straight test section, was of sufficient length to 
ensure that the flow was hydrodynamically developed for 
the flow velocities investigated. Based on these assump- 
tions, the experimental configuration was modeled as 
forced convection between parallel plates, for hydro- 
dynamically developed and thermally developing flows, 
subjected to a constant heat flux. The boundary conditions 
are one wall at a constant specified heat flux, and the 
other wall insulated. 

To analyze this condition, several quanties were 
defined as follows: 

The average heat transfer coefficient between the 
heated wall and the flow of air was defined by the equation 

Q = hA AT 

air 
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where, 'Q - ' was the heat convected to the air, ’ E - ’ 

was the average heat transfer coefficient, ' A p g' was the 
area of the heated Temsheet for the section being tested, 
and 'AT* was the difference between the average heated 
wall temperature (T ) and the average bulk temperature 
of the fluid (Tg) . The average bulk temperature of the 
fluid was defined as the arithmetic mean of the fluid 
inlet bulk temperature (Tj^) and the fluid outlet bulk 
temperature (Tgyp) . The delta temperature (AT) could then 
be expressed as: 



AT = T - T = T 

wo B wo 



T + T 
IN * O UT 

2 



The actual heat convected to the air (Q . was 

^air) 

calculated separately using the expression: 



Q - m c (T - T ) 

air p air OUT IN 



where ’c . ' was the specific heat of the air at constant 

pair ^ 

pressure and 'm' was the mass flow rate of the air. The 
mass flow rate, is further defined as the product of the 
volumetric flow rate of the air and the local flow density 
of the air as it transits the rotameter. The volumetric 
flow rate of the air was corrected to standard conditions 



for the rotameter and the local flow density was calculated 
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assuming "perfect gas" behavior at constant pressure and 
utilizing the temperature of the air as it exited the 
channel • 

The average Nusselt number was then calculated from the 
average heat transfer coefficient, with the expression: 



in which ’d* is the height of the channel and ’K . ' is the 

air 

thermal conductivity of the air evaluated at the average 
bulk temperature. 

The Reynolds number was then calculated for each test 
run as follows : 



where again *ra’ and ’d* are the mass flow rate and height 
of the channel respectively, 'A ' is the cross-sectional 
area of the channel, and ’y a p r ' is the dynamic viscosity 
of the air evaluated at the exit flow temperature. 

For the experiment runs utilizing the curved section, 
the Dean nuraher, defined as: 



h d 



Nu 




m d 



Re 



d 




De = Re 
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was also evaluated. In this equation, *d’ is again the 
height of the channel, and ' Rj_ r is the radius of curvature 
of the unheated wall surface of the channel (the convex 
wall) . 

The heat losses, by conduction, through the heated 
plate and the insulation, by conduction, were also cal- 
culated for each test section utilizing the expression: 

AT 

INS 

Q = 

lo ( x /K ) A pT 
INS INS 



With T AT j jsj g ' being the difference in temperature, as record- 
ed, between the first and second layers of AkMAPLEx 22 
insulation, and that recorded between the second and third 
layers. was ^he thickness °f each layer of insula- 

tion, 1 ^ xms * was t ^ ie thermal conductivity of the 
insulation, and ' Ap^* was again the area of the heated 
surface . 



In order to verify that Q • , as calculated by the 

3. ir 

temperature rise in the air, was actually the convection 
heat transfer from the Temsheet heated surface, an 
estimate of the radiation heat transfer was also calculated 
for each test run. The eight thermocouples in the lower 
Plexiglas surface measured the temperature of that surface 
and allowed the computation, utilizing the Stef an-Bolt zman 
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law of radiation, of the radiated heat transfer (Q ) with 
the expression: 



Qr 



o' (T 



4 

wo 



R 

R 



T 




where *cr ’ is the Stef en-Boltzman constant, and 1 R 1 is the 
total resistance to radiation heat transfer between the 
surfaces as computed by: 



R 



R 





'wo 



e . 

Wl 



1 ) 



In the above equation, ' e 'is the emissivity of the 
n ’ wo 7 

heated Temsheet wall surface and 'e . ' is the emissivity of 

wi 7 

the unheated adiabatic wall surface (Plexiglas) , 

A sketch of the control volume and the energy balance, 
with sample calculations for one test run, is given in 
Appendix B. 



B. RESULTS 

The data obtained from the experimental test runs were 
evaluated utilizing the expressions described in the 
previous section, and as illustrated in Appendix B. The 
results of the major variables are shown in tabular form, 
in Table I for the straight test section, and Table II for 
the curved test section. A plot of the average Nusselt 



number verses Reynolds number, is given in Figure 9 for 
comparison of the two test sections. The uncertainty bands 
have been indicated for the plotted experimental data and 
a sample calculation for the complete error analysis is 
given in Appendix A. 

As shown, the results indicate an increase in the 
rate of heat transfer with increasing Reynolds number for 
both the straight and curved test sections. Additionally, 
the heat transfer rate of the curved test section was 
higher than that of the straight test section for each 
Reynolds number investigated. 

In investigating the results obtained, and considering 
those Reynolds numbers below 2000 (1025 based on this 
channel configuration) as laminar, at least squares 
straight line correlation of the data results in: 

0.67 

Ru = 0.06 5 Re d 

for laminar flows, and 

— 0.58 

Nu = 0.117 Re d 

for transition regime flows, in the curved test section. 
With correlations of 



0.25 

Nu = 0.983 Re 
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for laminar flows, and 



— 0.56 

Nu = 0.110 Re d 

for the transition flows in the straight test section. 

In the laminar flow regime, the increase in heat 

transfer was approximately fifteen percent, while in the 

higher transition flow region the increase was approximately 

thirty percent. This compares favorably to the earlier 

studies by Kreith .[Ref. 4], in which he reported that the 

heat transfer rate along concave walls increased from 

twenty-five percent to sixty percent for Reynolds numbers, 

4 6 

based on hydraulic diameter, between 10 and 10 , and 
Ballard [Ref. 30] , who found an eleven percent increase 
in heat transfer rate for laminar flows. 
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Reynolds Number 

Figure 9. Straight vs. Curved Section results for the present Study 



Re d 

784 

941 

1099 

1257 

1414 

1571 

1729 

1886 

2042 

2203 

2361 

2517 

2675 



TABLE I 



SUMMARY OF STRAIGHT TEST SECTION RESULTS 



Qair FT AT Nu 

(W) (W/m 2 ° C) (°C) 



39.42 


21.28 


40.14 


22 . 56 


43.26 


23.11 


45.13 


23.93 


51.80 


25.89 


50.44 


26.59 


53.50 


28.61 


57.01 


31.32 


59.24 


31.99 


61.65 


33.50 


62.97 


35 . 30 


67.11 


35.44 


69.14 


36.39 



26.7 


5.08 


24.0 


5.38 


25.2 


5.51 


25.4 


5.71 


27.0 


6.18 


25.6 


6.34 


25.2 


6.83 


24.5 


7.47 


25.0 


7.63 


24.8 


7.99 


24.0 


8.42 


25.5 


8.45 


25.6 


8.68 
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Re d 

782 

937 

1092 

1172 

1251 

1406 

1570 

1720 

1883 

2038 

2196 

2351 

2505 

2659 



TABLE II 



SUMMARY OF CURVED TEST SECTION RESULTS 



De. 


Qair 

(W) 


h 

(W/m 2 °C) 


f — \ 

H U 
<3 o 

V — f 


Nu 


114 


41.19 


23.63 


24.3 


5.64 


137 


47.18 


26 . 57 


24 . 7 


6.34 


159 


50.88 


28.56 


24.8 


6.81 


171 


53.71 


30.49 


24.5 


7.27 


183 


57.16 


32.32 


24.6 


7.71 


205 


58.11 


33.01 


24.5 


7.88 


229 


6 2.63 


35.12 


24.9 


8.38 


251 


69.41 


37.23 


26.0 


8.88 


275 


70.90 


38.29 


25.8 


9.13 


297 


76.55 


41.66 


25.6 


9.94 


320 


80 . 36 


42.72 


26.2 


10.19 


343 


83.30 


45.98 


25.2 


10.97 


366 


88.82 


46.87 


26.4 


11. 18 


388 


94.67 


49.88 


26.4 


11.90 
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V. DISCUSSION AND CONCLUSIONS 



With. the measurement of the unheated wall and the air flow 
temperatures, it was indicated that there was a negligible 
difference between the fluid bulk temperature and the tempera- 
ture of the unheated wall. Additionally, the energy balance 
and sample calculations, provided in Appendix B, demonstrate 
that the radiated heat transfer from the unheated wall sur- 
face is minimal. Thus it was assumed that the heat transfer 
to the flowing air was solely by convection from the heated 
Temsheet . 

The high aspect ratio of the channel as liscussed pre- 
viously, provided the original basis for the assumption that 
the experimental apparatus, as configured, was one of infinite 
parallel plates. The experimental data, for those Reynolds 
numbers in the laminar flow region, tend to substantiate this 
assumption, in that they approach the theoretical limit for 
average Nusselt number of 5.385, for parallel plates with one 
wall heated at a constant heat flux and the opposite wall 
adiabatic [Ref. 37]. For a Reynolds number, based on hydraulic 
diameter, of 1835, the average value of Nusselt number obtained 
was 5.382. 

With the formulation of the problem complete, the 
experimental data obtained can now be compared to the analytical 
solutions and experimental results of the same or similar 
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problem. That problem being, flow between infinite parallel 
plates with one wall at constant heat flux and the opposite 
wall adiabatic. 

A comparison of the experimental results obtained in this 
study and the experimental results of Durao [Ref. 29] and 
Ballard [Ref. 30] (who investigated the problem solely for 
laminar flows), are shown in Figure 10, for the straight sec- 
tion, and Figure 11, for the curved section. In both these 
studies, as in this investigation, the channel has an aspect 
ratio of 40, and the experimental procedures were similar. 

It can be seen that the results obtained in this study, in 
the laminar flow regime, appear to be in agreement with the 
data obtained previously for both test sections. Additionally, 
numerical correlations of the present data, as given earlier, 
for the laminar flows compares favorably to the numerical 
correlations of both Durao and Ballard, for Reynolds numbers 
greater than 600, as given in Appendix C. In the two investi- 
gations, mentioned above, comparisons were made to the analy- 
tical studies of McCuen, et. al. [Ref. 35], for heat transfer 
between infinite parallel planes with constant wall temperatures 
and heat flux. Ballard also compared his data to the 
Worse-Schimdt solution [Ref. 36] , for heat transfer in tubes 
and annular passages. In both the Durao and Ballard studies, 
the experimental data compared relatively well to the predicted 
solutions of MCuen, et al. A curve of the predicted solution 
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of the McCuen, et. al . study is shown in Figure IQ. The re- 
strictions of that study were as follows: 

1. the velocity profile is fully established. 

2. the fluid transport properties and density are assumed 
constant . 

3. axial heat conduction is negligible. 

4. viscous energy dissipated is negligible. 

5. conditions are invarient with time (steady state). 

All these restrictions have been accounted for in the 

design of the experimental configuration and/or verified in 
the experimental processes. A dimensionless axial length 
coordinate defined as: 

x * = * 

D^Re Pr 

was calculated for 'x' equal to 29.2 centimeters (the length 
of the straight test section). In this equation, 'Pr ' is 
the Prandtl number of air, taken as 0.7, ' ' is the hydraulic 

diameter Cl .239 centimeters), and 'Re ' is the Reynolds number, 
based on the hydraulic diameter. 

Based on the values of x* , it was determined that the 
experiments, for laminar flows, were in the thermal entrance 
region. It was therefore decided to integrate the expression 
of the local Nusselt number as given by McCuen, et. al. [Ref. 
35] : 

0.67095 x*^-> - 2 Cl +^i)x* 

Q 

wo 
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for the small values of x, in order to obtain the following 
expression for the average Nusselt number: 



Nu = - 



4(1 * ^i)x* 
^wo 



ln (1 ‘ 0.67095 + Q W1 ^ X * 

x wo 



Q,.,i -( 2 / 3 ) 



) 



In this expression, Q • was the heat convected to the fluid 

W1 

from the unheated wall, and Q was the heat convected to the 
fluid from the heated wall. In accordance with the previous 
assumptions, was considered equal to zero, and Q^. q was 

considered equal to Q a ^ r > an d the average Nusselt number could 
then be calculated, and the expression reduces to: 



rj — -0.75 , f , 

Nu = —£ — In (1 - 



.*-( 2 / 3 )) 



x- — 0.67095 

This value of average Nusselt number was plotted for values 
corresponding to an x* less than 0.042, as this represents the 
maximum value for which the analytic solution for Nusselt 
number is valid. This value of x*, (0.042), for x equal to 
29.2 centimeters, equates to a Reynolds number value of 
approximately 800, based on hydraulic diameter (or 411 based 
on the channel height) . 

The graphical representation shows that the experimental 
results of this study as well as those of Ballard and Durao 
plotted above the analytical solution. This difference can 
probably be explained, at least in part, by the difference in 
geometries used in the studies, and by the limitations inherent 
in any experimental work. The side wall effects and the in- 
ability to totally account for all the heat transfer processes 
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Figure 10. Comparison of present data w i tli Durao, Ballard and McCuen for Laminar 
Flows, Straight Section 
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Figure 11. Comparison of present data with Durao and Ballard for Laminar Flow, 
Curved Section 



and/or losses could account for the differences of approximately 
nine to ten percent. 

In the comparison of the experimental results obtained in 
this study for those flows not in the laminar regime, the 
comparison and correlation was made to those analytical solu- 
tions and experimental results for completely turbulent flows. 
For the most part, the analytical solutions are valid for 
Reynolds number values greater than 10 4 (equating to a Reynolds 
number of 5123 in this study) , and tube flow (or other geo- 
metries based on hydraulic diameter) . 

A graphical comparison of the results obtained in the 
straight test section runs with the analytical solutions and 
other experimental results is made in Figure 12. The plot 
shows the curve of heat transfer in a straight tube for con- 
stant wall temperature, as given by the Dittus -Boelter equations 
[Ref. 38]: 

Nu = 0.023 Re 0,8 Pr 0,4 

and by the empirical relation as given by E.R.G Eckert [Ref. 

39] : 



Nu » 



0.0396 Re °' /5 Pr 



(1 - 



Re 



44 , 

0 . 1 2 5 J 



This analytical expressions when extrapolated down show some 
agreement with the values obtained at flow rates between forty 
and seventy percent of the rotameter full scale. Unfortunately, 
those values obtained at the higher flow rates showed less of 
a correlation than those mentioned above. These results, 
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obtained for volumetric flow rates higher than seventy percent 
may have been affected by the fact that they were near the top 
end of the rotameter scale, and therefore, may inherently be 
somehwat less accurate than the mid-scale readings. The ther- 
mocouple problem discussed previously, may also have contri- 
buted to the decrease in correlation. 

When extrapolated down, the experimental results of Kays 
and Leung [Ref. 28], for heat transfer in annular passages, 
with one wall heated at a constant heat rate, and r* equal to 
1.0, (where r* is the annulus radius ratio r^/r Q ) , equivalent 
to parallel plates, appears to correlate well with the present 
results. This extrapolation, was again utilizing a least 
squares straight line correlation, as given in Appendix C. 

Also shown are the analytical solutions of Petuhkov and Popov 
[Ref. 25] , for turbulent heat transfer in round tubes given 
as : 



where * f * the friction factor, has been calculated using the 
Filonenko equation given in Karlekar and Desmond [Ref. 40] as: 
f = (1.82 log Re - 1.64) ' 2 

and , 

k x (f) = 1 + 1 . 4f 



Nu 



(f / 8) Re Pr 



k 1 (f) + k 2 (Pr) /f7"8 (Pr^ 2 ^ 



1 ) 



k 2 (Pr) 
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The analytical expression, given by Shibani and Ozisik [Ref. 
26], for turbulent heat transfer between parallel plates: 

Nu. = 12 + 0.03 Re a Pr b 



where , 

a = 0.88 - 



0.24 

(3.6 + Pr) 

-0.6 Pr 



b = 0.33 + 0.5 e 
valid for 0.1<Pr <10, and 10^<Re <10 b , is also plotted. 

While no conclusions can be drawn from the data obtained 
because it is not in the range of Reynolds numbers of these 
studies, it is worthy of note, that the present data as extra- 
polated, and the analytical and experimental works of others 
can be graphically shown to be close. It should also be noted 
that the investigations of Hsu [Ref. 34], Ashton [Ref. 33, and 
Haynes and Ashton [Ref. 32] all plotted experimental results 
between the Petuhkov and Popov expression and that of Shibani 
and Ozisik, as shown graphically in the Haynes and Ashton 
study above. 

The curved section data is plotted in Figure 13, and 
compared to the experimental work of Brinich and Graham [Ref. 
31] for flow and mass transfer in a curved channel. The 
analytical solution of Petukhov and Popov, and Shibani and 
Ozisik are also plotted as a reference. A least squares 
correlation as given in Appendix C was used to extrapolate 
the data of Brinich and Graham. The accuracy of the data 
points used in this correlation are subject to errors, in that 
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the actual values were not given in their study and had to be 
taken from a plot of Stanton number versus Reynolds number as 
given in the reference. 

Irrespective of the relatively small differences in the 
results and the correlations, the trends appear to be conclu- 
sive. The heat transfer rates can be seen to increase with 
increasing Reynolds number, and also that the heat transfer 
rate in a curved channel is measurably increased over that 
rate in a curved channel is measurably increased over that 
rate in a straight channel. This increase in heat transfer 
rate, appears to begin with the development of the Taylor- 
Goertler vortices and continues into the transition and tur- 
bulent flow regions. This result also appears to have been 
reached by K. V. Dement ’ eva and I. Z. Aronov [Ref. 41], who 
in 1978 , investigating hydraulic resistance and heat transfer 
in curvilinear channels of rectangular cross-section, 
concluded that there was an enhancement of heat transfer 
intensity in curved rectangular corss - sect ion channels as 
compared to straight pipes and channels. 
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Reynolds Number 

Figure 13. Comparison of Present Data for Transition Flows, Curved Section 



VI. RECOMMENDATIONS 



The amount of experimental work remaining in the study of 
heat transfer in curved channels is significant. Additional 
experiments in the transition region and in the fully turbu- 
lent region of fluid flow, should be conducted to refine the 
correlations for average Nusselt number. With this work, it 
would be necessary to design and construct a larger experi- 
mental configuration, capable of higher volumetric flow rate 
and greater temperature differences between the heated surface 
and the fluid inlet temperature. Additionally, care should 
be taken to ensure that the temperature of the exit flow, of 
each test section, is accurately measured. This may be 
accomplished by the use of a mixing cup arrangement at the 
exit flow locations. In a channel designed and constructed 
as the one in this study, care must be taken to allow for the 
removal of the mixing cup at the straight section exit location, 
to ensure that turbulence is not induced in the flow by the 
thermocouples prior to the fluid entering the curved section. 

Also in the way of temperature measurement, a thermocouple 
should be installed at the entrance to the rotameter, or pre- 
ferably some better device or means of determining the 
volumetric flow rate. This would ensure a more accurate 
measurement of the fluid temperature for an evaluation of the 
fluid properties in the mass flow calculation. 
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APPENDIX A: EXPERIMENTAL UNCERTAINTY 



The uncertainties for the major variable in the experiments, 
were calculated in accordance by the method described by S. Kline 
and F. McClintoch [Ref. 42] . The estimates of the uncertainty 
in the measured quantities were made conservatively, and as such, 
there is considerable confidence in the uncertainties as calcu- 
lated. An example of the calculations follows. 

The Reynolds number was defined by the equation: 

Ke d p . A 
M air c 

As such, the uncertainty was calculated by : 
dRe 



Ri" 



- 

▼ m 



<v * (^) 2 * (^i £) 2 * qV- 

m a y air A c 



The uncertainty in the mass flow rate was determined in a 
similar manner, by its equations: 



m = 



and 



- Q I 

RT 



4 - { 

m 



c £) 2 * #) 2 * #) 2 - C & 2 



The uncertainties in the height of the channel, dynamics vis- 
cosity of air, cross-sectional area of the channel, the gas con- 
stant, pressure, temperature and the reading of the volumetric 
flow rate were obtained from estimates as: 



dd 

i ~ 


.01987, 


dp„ . 
air 

\i . 


dA c 

= .00223, , c 

A 


= .01976 






air 


c 




pu 

|*j 

ii 


.00002, 


dl pu 

ii 


.00048, = 


.00055 
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= .0200 

The uncertainty in the mass flow rate was then calculated as; 
dm 



= .0200 



m 



and the Reynolds number uncertainty was 



dRe ( 

ReT 



- Vc. 



0 200) 2 + (. 01987) 2 + (.0022. 3 ) 2 + (.01976) 2 = .0345 



Re d = 2042 + 71 

The uncertainty in the Nusselt number was also calculated 



in an identical manner, using the equation: 

rj — Ed 

Nu = 



K . 
air 

and the uncertainty as 



dNu 

Nu 






_■) * c^r * c^ii ) 2 

h d K air 



It was therefore, necessary to first calculate the uncertainty 

for h. Using the same approach, 

E = Qair 



A PL iT 



and 



dh _ jA Qa 

h " V^ ai 



air, 2 , r dA PL, 2 t ,dAT, 2 
air J L A T , t J l AT J 



TL 



where AT = T 



wo 



B 



dA 



A 



— = .0026 and = .0048 



PL 



AT 



The uncertainty in Qair was calculated by 



Qair m C pair (T qut - T JN ) 
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and 



dQa 



Qair 
where 



air _ J r dm^2 , , dC Pair, 2 , r d(T OUT T IN),2 

— - t/ (.— rJ + l-p J + (.t r - T J 

ir 1 m ^Pair OUT IN 



— = 0.0200 , 

• 7 

m 



dC 



Pair 



'Pair 



= .0042, 



and d ( T QUT ~ t in ) 
T OUT " t in 



= .0623 



These calculations result in the following uncertainties: 
- .0656 

Qair 

— = .0658 



h 

dNu 

Nu 



= .0688 



It is therefore, the uncertainty in the difference between 
outlet and inlet temperatures of the flowing air, that is 
ultimately the major source of average Nusselt number uncer- 
tainty. This then is the basis for the recommendation to 
better ascertain the exit temperatures for each test section. 
The calculations, as shown above, were for the 65% straight 
section run. Similar calculations for the curved section at 
65% yield the following; 
dRe^ 

= .0345 



Re 



Re d = 2038 + 70 

pE = .0563 
Nu 
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Nu = 9.34 +_ 0.525. 

The values for the uncertainties for all the variables in 
the curved section test run are as follows: 



QUANTITY 


UNCERTAINTY 


a pl 


.0026 


Ac 


.0198 


c p 


.0042 


d 


.0199 


De 


.0398 


h 


. 0526 


Kair 


.0003 


Nu 


.0563 


• 

m 


. 0200 


Q 


. 0200 


Qair 


.0523 


Re 


. 0345 


t b 


.0073 


t in 


. 0029 


t out 


.0106 


T w 


.0060 


T -T 

1 OUT IN 


.0482 


AT 


.0047 


U 


.0022 


P 


.0007 
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APPENDIX B: SAMPLE CALCULATIONS 



A sketch of the control volume for the energy balance, 
for the straight section, indicating the major heat transfer 
components is shown in Figure 14. The sample calculations 
that follow, show the methods and equations that were used 
to obtain these heat transfer components. The heat transfer 
components were required to satisfy the energy balance, and 
to obtain an average Nusselt number as a function Reynolds 
number of the flowing air. 

The curved section energy balance and calculations are 
similar 




Temsheet 



FIGURE 14. Energy balance in straight section 
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A. SAMPLE CALCULATIONS DATA 



T ROOM 

Q 

V PR 

V H 

r pr 

t in 

t insi 

T INS2 

Twol = 

Two 2 = 

Two 3 = 

Two4 

Two5 = 

Two6 = 

Two 7 = 

Two 8 = 

Twil 

Twi2 

Twi3 

Twi4 = 

Twi5 

Twi6 

Twi7 

Twi8 = 

T 0UT 

t exit 



21 


. 6° i 


n 


55 






3. 


216 


5 V 


37 


.1210 V 


2. 


078 




0 . 


839 


mV 


1 . 


324 


mV 


1 . 


583 


mV 


1 . 


840 


mV 


1 . 


840 


mV 


1 . 


987 


mV 


2. 


001 


mV 


2. 


112 


mV 


2. 


125 


mV 


2 . 


245 


mV 


2. 


241 


mV 


0 . 


964 


mV 


0 . 


967 


mV 


1 . 


012 


mV 


1 . 


008 


mV 


1 . 


056 


mV 


1 . 


047 


mV 


1 . 


101 


mV 


1 . 


100 


mV 


1 . 


124 


mV 


1 . 


076 


mV 


0 . 


074 


^ 2 
2 m 



21 .82°C 
33 .34°C 
39 . 45° C 



Two = 



Twi = 



= 2 7 . 9 7° C 

= 26 . 91°C 



50 .09°C 



25. 76°C 



= 80 . 44 0 F 
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